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Abstract 



"t^ ' We define a deterministic "scattering" model for heat conduction which is contin- 

uous in space, and which has a Boltzmann type flavor, obtained by a closure based 
on memory loss between collisions. We prove that this model has, for stochastic 
driving forces at the boundary, close to MaxwelUans, a unique non-equilibrium 
steady state. 

^ I 1 Introduction 

O 

cn I In this paper, we consider the problem of heat conduction in a model which is a 

^ ■ continuous approximation of a discrete model of a chain of cells, each of which 

contains a (very simple) scatterer in its interior. Particles move between the cells, 

Q I interacting with the scatterers, but not among themselves, similar to the model put 

forward in [3]. 

After a quite detailed description of this scatterer-model, we will finally arrive 
in Eq. (3.5) at a formulation with a continuous space variable x varying in [0, 1]. 
C^ I This approximation will be obtained by taking formally the limit of A^ ^ oo cells, 

but taking each cell of length 1/A^. (The reader who is interested only in the formu- 
lation of the a:;-continuous equation can directly skip to Eq. (3.5).) We then proceed 
to show our main result, namely the existence of solutions to this Boltzmann-like 
equation for initial conditions close to equilibrium. Therefore, this will show that, 
although the model is deterministic (except for the boundary conditions), with no 
internal dissipation, it has a (unique) non-equilibrium steady state when driven 
weakly out of equilibrium. The only approximations of the model are the limit of 
N —>■ oo, and a closure relation which models a loss of memory between collisions. 

1.1 One cell 

To define the model, we begin by describing the scattering process in one cell. We 
begin with the description of one "cell". The cell is 1 -dimensional, of length 2L, 
and with particles entering on either side. These particles have all mass m, velocity 
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V and momentum p = mv. These particles do not interact among themselves. 
Note that v G M and more precisely, w > if the particle enters from the left, 
while f < if it enters on the right side of the cell. In the center of the cell, we 
imagine a "scatterer" which is a point-like particle which can exchange energy and 
momentum with the particles, but does not change its own position. (This scatterer 
is to be thought of as a 1-dimensional variant of the rotating disks used in [3].) The 
scatterer has mass M and its "velocity" will be denoted by V. The collision rules 
are those of an elastic collision, where v and V denote quantities after the collision 
while V, V are those before the collision. In equations, 

i = -Qv + (1 + q)V , 
V ={l- q)v + qV , 

with 

M — m m 1 — f? 

M + m M 1 + Q 

Note that g G (—1, 1), since we assume m and M to be finite and non-zero. If 
i; > 0, we say that the particle leaves the cell to the right; if i; < 0, we say it leaves 
to the left. 

For simplicity, we will assume ^ > 0, that is, M > m. For the momenta, we 
get the analogous rules 

P = -QP + C^- Q)P , 
P ={1 + q)p+qP. 



Note that the matrix 



« = ir/, Ti <'-^' 



has determinant equal to —1 and furthermore S'^ = I. 

We next formulate scattering in terms of probability densities (for momenta) 
for just one cell. We denote by g{t, P) the probability density that at time t the 
scatterer has momentum P (= MV) and we will establish the equation for the time 
evolution of this function. To begin with, we assume that particles enter only from 
the left of the cell, with momentum distribution (in a neighboring cell or a bath) 
p 1-^ /]^(i, p), where p = mv. Thus, there are, on average, pf^{t, p)Ap/m particles 
entering the cell (per unit of time) from the left with momentum in \p,p + dp\. 
Note that f^ has support on p > only, (indicated by the exponent "+"); it is 
the distribution of particles going to enter the cell from the left. Also note that the 
distribution of the momenta after collision, i.e. , before leaving the cell, is in general 
not the same as f^. 

Denote by ^(t) the stochastic process describing the momentum of the scat- 
terer. We have for any interval (measurable set) of momenta A, for the probabilities 
P: 

P{^{t + dt) £ A) = P(^(t) £ A; no collision in [t, t + dt]) 
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+P(^(t + dt) £ A; collisions occurred m[t,t + dt]) . 

We assume for simplicity that with probability one, only one collision can occur in 
an interval [t,t + dt]. If there is a collision in [t, t + dt] with a particle of velocity 
V = p/m > 0, this particle must have left the boundary at time t — ^L with 
momentum p. Therefore, 

P{^{t) £ A ; a collision occurred in [t,t + dt]) 

= dt [ dP [dp5{P -gP-{l + Q)p) g{t, P) y+{t - fL,p) . 

This equation neglects memory effects coming from the fact that a particle may 
have hit the scatterer, bounce out of the cell and reenter to hit again the scatterer. 
Similarly, 

P(=^(t) G ^ ; no collision occurred in [t, t + dt]) 

= (^l-dtljp^f+{t-fL,p)^ J^dPg{t,P). 
We immediately deduce the evolution equation, 

dtg(t, P) = - g{t, P) [dp ^f+{t - fL,p) 

Jr+ 



+ ^ldpgit,^^^^^)^fMt-^L,p) 



(1.3) 



Q.'-- ' 



Note that this equation preserves the integral of g over P, i.e., it preserves proba- 
biUty. 

This identity generalizes immediately to the inclusion of injection from the 
right, with distribution f^ having support in p < 0. One gets 



(1.4) 



dtg{t,P) = -g{t,P) / dp^[f+{t-fL,p) + f^{t + fL,p) 

JR ^ 

+ i / dpgit, ^^M^) M (/+(i _ fL,p) + f-it + fL,p) 
Q Jr i; \ t- 

In the stationary case, this leads to 

g{P) = - [ dpg{^^^^) M (/+(p) + f-(p)) , (1.5) 

Q^ JR 



where 



A= /dpM (/+(p) + /^(p)) (1.6) 

JR 

is the particle flux (see Sect. 3.1 below). 

It is important to note that the solution g of Eq. (1.5) only depends on the sum: 
f = f]^ + f^, and thus, we can define a map 
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where gf is the (unique) solution of Eq. (1.5). It will be discussed in detail in 
Sect. 5. 

We can also compute the distribution of the momenta of the particles after 
collision. We have 

P(p G ^ ; a collision occurred in [t,t + dt]) 

= dt f dp f dP6{p + Qp-{l-g)P)g{t,P)^-^f{t-^L,p). 

This particle reaches the left or right boundary of the cell (according to the sign of 
p) after a time mL/\p\ (assuming the scatterer is located in the center of the cell). 
Therefore, we have for the ejection distributions /^J (on the left) and f^ (on the 
right): 

S/i:(t,P) = ^ /dp<,(t- ^L,^)M;(t_ ™L- ^L,p) , 
and 

Mim =^-^ Ijpgit- ^L/E^)\^f(t- ^L- ^L,p) , 
where 9 is the Heaviside function. In the stationary case we get 

S/i:(P) = ^/dP5(^)S/(^), (1.7) 



and 



S/r^® = f^ / dp <7(^) S/(P) . (1.8) 



Since g = gj is determined by the incoming distribution /in = /l^ + /r (and is 
unique if we normalize the integral of ^f to 1) 



/ 

JR 



dPg{t,P) = l, (1.9) 



we see that the outgoing distribution /out = /l + /r^ is entirely determined by the 
incoming distribution. Note also that the flux is preserved: 



y"dpg/i„(p) = y"dpM/out(p). 



1.2 Stationary solutions for one cell 

Here, we will look for stationary states of the evolution equation (1.3), which have 
also the property that the ejected distributions are equal to the injected ones. It is 
almost obvious that Maxwellian fixed points can be found, but for completeness, 
we write down the formulas. The reader should note that the distributions /in and 
/out have singularities at p = 0. This reflects the well-known fact that slow particles 
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need very more time to leave the cell than fast ones. However F(j)) = —f{p) is a 
very nice function, and it is this function which appears in all the calculations of 
the fluxes, and stationary profiles. In this section we do the calculations with the 
quantity /. Starting from Sect. 3, we will use F. 
We impose the two incoming distributions 

/l+(p) = ^e(+p)^e-'3pV(2m) ^ 
and 

where a is an arbitrary positive constant (related to A in (1.6)) and 6 is the Heaviside 
function. It is easy to verify, using Gaussian integration and the identity M = 
Mq^ + m(l + of, that the solution of equation (1.5) is given by 



/5 „-/3p2/(2M) ^ / 1^ -/3p2(l-p)/((l+p)2m) 



/+(p) = ^e(+p) ™e-^PV(2m) _ 



Moreover, using the same identity several times one gets from Eqs.(1.7) and (1.8) 
for the exiting distributions 

/iTCp) = a0(-p)^e-^PV(2-) ^ 
and 

/R \P) — "l^V-T/^;]^ 

Therefore, we see that the Maxwellian fixed points (divided by \p\) preserve both 
the distribution g of the scatterer, as well as the distributions of the particles. 
In fact, there are also non-Maxwellian fixed points of the form 

/+(p) = a0(+p)^e-^(^-™«)'/(2™) ^ 

and 

/r (P) = a0(-p)^e-'3(^-™«)'/(2™) 

It is easy to verify that the solution of equation (1.5) is now given by 



f^ „-/3(P-Ma)2/(2M) 



Moreover, 



and 



^^^^ = ^/2^^ 



/l-(p) = a0(-p)^e^'5(^-™")'/(2™) 



/+(j9) = a0(+p)f^e-'5(P-™")'/(2m) 



The verification that this is a solution for any a € M is again by Gaussian integra- 
tion. Note that if a / there is in fact a flux through the cell. 
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2 A^ cells 

The model generalizes immediately to the case of N cells which are arranged in a 
row, by requiring that the exit distributions of any given cell are equal to the entry 
distributions of the neighboring cells: The cells are numbered from 1 to A^ and we 
have the collections of functions f^-, f^^, f^^, and f^- for the particle fluxes and 
gi for the scatterers, i = I, . . . ,N. The equality of entrance and exit distributions is 
given by the identities f^^^i = f^^, and /j^^ = /^Jj^^ for 1 < i < A^. The system 
is completely determined by the two functions f^^ and f^^- The equations (1.5) 
generalize to 



g.{P) = ^ / dp5.(^^^^™)f ( ft.ip) + f^P) ) , (2.1) 



.P-(1+£)pn IpI ( f+ 



and similarly (1.7) and (1.8) lead to 



>^ — Q JM. 
>^ — Q JM. 



where fi = /£j"^ + /^^. Clearly, the Gaussians of the previous section are still 
solutions to the full equations for N contiguous cells. 

Here we have closed the model by assuming independence between the parti- 
cles leaving and entering from the left (and from the right). In concrete systems 
this is not true since a particle can leave a cell to the left and re-bounce back into 
the original cell after just one collision with the scatterer in the neighboring cell, 
and, in such a situation there is too much memory to allow for full independence. It 
is possible to imagine several experimental arrangements for which independence 
is a very good approximation, see also [3, 2] for discussions of such issues. One of 
them could be to imagine long channels between the scatterers where time decor- 
relation would produce independence. Note that "chaotic" channels may be more 
complicated since they can modify the distribution of left (right) traveling particles 
between two cells. 

3 Continuous space 

We are now ready to formulate the model in its final form. The cells are now 
replaced by a continuum, with a variable x G [0, 1] and the relations we have de- 
rived so far will be generalized to this continuum formulation. So we have moving 
particles, of mass m and described by a time-dependent density f{t, p, x). 

The scatterers have mass M and their momentum distribution is called g{t, P, x). 
It is best to think that the continuous variable x G [0, 1] replaces the discrete index 
i G {0, . . . , N}. There is then an implicit rescaling of the form x « i/N. Recall 
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that the scatterers are fixed in space (ahhough they have momentum) but that the 
particles will move in the domain [0, 1]. 

We first impose, for all x G [0, 1], the normalization 



dPg{t,P,x) = l, (3.1) 

which is the generalization of Eq.(1.9). The particles again do not interact with 
each other, but only with the scatterers and, expressed in momenta, the matrix 
maps (p, P) to (jj, P): 

as in Eq.(1.2). 

We now rewrite the problem in the form of a Boltzmann equation, which takes 
into account this matrix, as well as the particle transport. One obtains, with (p, P) 
related to {p, P) as above: 

dtfit,p,x)-\ dj:f(t,p,x) 

m 

dP (^-^f{t,p,x)g{t,P,x)-^-^f{t,p,x)g{t,P,x)) 
dtg{t,P,x)= [dp (^-^f{t,p,x)g{t,P,x)-^-^f{t,p,x)g{t,P,x) 

The time independent version of the equation will be derived below from the model 
with a chain of cells. It is useful to introduce the function 

F{t,p,x)=^^f{t,p,x), 
m 

and then Eq.(3.3) takes the form 

mdtF{t, p, x) + pdxF{t, p, x) 

dP (F{t, p, X) g{t, P, X) - F{t, p, X) g{t, P, x) 



(3.4) 
dtg{t, P,x)=J dp (F{t,p, X) g{t,P,x)- F{t,p, x) g{t, P, x)) . 

Remark. One can also imitate a scattering cross section by introducing a factor 
7 G [0, 1] in Eq.(3.4) (in the integrals) but this can be scaled away by a change of 
time and space scales. See also Sect. 8.3. 

We come now to the main equations whose solutions will be discussed in 
detail in the remainder of the paper. The equation (3.4) takes, for the stationary 
solution, the form 
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0= dp (^F(p, X) g{P, X) - F(p, X) g{P, x)J . (3.5b) 

We will show that this equation has non-equilibrium solutions. 
Remark. Note that the model we have obtained here is not momentum-translation 
invariant, because of the term sign(j)), except when the r.h.s. of the equation is 0. 
Derivation of (3.5). The derivation of (3.5) from (2.1-2.2) is based on the follow- 
ing formal limit: We replace the index i by the continuous variable x = i/N and 
set e = 1/{2N). We consider that /^- is at {i - \)/N = x - e, while /^. is at 
x + e. We have the correspondences, with 0±(p) = 9{^p): 

d+(p)F(p, x-e) = M/+^ , e^(p)F(p, x-e)= M/^^. , 

e+(j>)F(p, x + e)= ^-^f+. , e-(j))Fip, x + £)= M/r;, , 

g{P,x) = g^{P). 
To simplify momentarily the notation, let 

F_(p, x) 9+{p) = e+(p)F{p, X - e) , F^ip, x) 9-{p)= 9^{p)F(p, x - e) , 
F+{p, X) 9+{p) = 9+{p)F{p, x + e), F+{p, x) 9-{p) = 9^{p)F{p, x + e) . 

With these conventions, (2.1) becomes (setting A = 1): 

g{P,x) = - I dqg{^^^^^,x)iF.{q,x)9+{q) + F+{q,x)9^{q)) , (3.6) 
which is equivalent to (3.5b). Similarly, Eq.(2.2) leads to 

F.(jj,x)9.(p) = ^^ / dqg{S±^,x){F.{q,x)9+{q) + F+{q,x)9-{q)) , 

0(-\-n) r (3.7) 

F+(j>,x)9+(j;>)=-^ / dqgii^,x){F.{q,x)9+{q) + F+{q,x)9.{q)). 

^ — Q Jk 

Subtracting the first equation from the second in (3.7) leads to 
F+(p,x)9+{p) - F^{p,x)9^{p) 

= '-^ f dqg{S^,x) {F.{q,x)9+{q) + F+{q,x)9.{q)) . ^^'^^ 
^- Q Jr ^ 

On the other hand, since dpdP = dpdP we can, by (3.5b), impose the condition 

dPg{P,x) = l, 
for all X. Then we have the trivial identity 

F-(j),x)9+(j))-F+(j),x)9-(j)) = [ dP g{P,x){F.(j),x)9+(j))-F+(j),x)9-(j))) . 

JR 

(3.9) 
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Subtracting (3.9) from (3.8), we get for p > 0, 



F+(p,x)-F_(p,x) =-^ f dqg{^,x){F^(q,x)e+(q) + F+(q,x)eAq)) 

^ ■' (3.10) 

while for p < 0, 

F+{p,x)-F^{p,x) = --^ fdqg(E±£^,x){F^iq,x)e+(q) + F+(q,x)e-iq)) 

y (3.11) 

Note now that 

F_{q,x)e+{q)+F+{q,x)e^{q) 

= F{q, X - e)9+{q) + F{q, x + e)e^{q) 
= F{q, X) + e+{q){F{q, x - e) - F{q, x)) 
+ e-iq){Fiq,x + e)-Fiq,x)). 

Therefore, replacing the F± in the r.h.s. in (3.10) and (3.11) by F{p, x) is a higher 
order correction in e, and we finally get 

F(p, x + e)-F(p,x-e) = ^-^^ j dq(g{P+M^ x)F{q, x) - g{q, x)F(p, x) 

A further change of integration variables leads to (3.5a), while (3.6) leads to (3.5b). 
(We have not taken into account the scaling by e = 1/(2 A'^) which is needed to get 
the derivative; we will come back to this question in the discussion in Sect. 8.3.) 
This ends the derivation of (3.5). 

The derivative term in Eq.(3.5) reflects the gradients which have to appear when 
the system is out of equilibrium. However, if the system is at equilibrium, the 
equivalence between Eq.(3.5) and Eqs.(1.5)-(1.8) immediately tells us that sta- 
tionary solutions in the form of Gaussians (for F, not for /) exist: 



P „-/3pV(2m) „(p^_ J P ^-0Py(2M) 



Furthermore, we have again translated versions of this fixed point, 



P „-/3(p-ma)2/(2m) „r r>-. _ J 1^ -0(P-Maf /{2M) 



0A3) 
because in this case, the r.h.s. of Eq.(3.5) is zero. 
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3.1 Flux 

We can define various fluxes of the particles (recall that F{p, x) = \p\f(jJ, x)/m): 



$P = particle flux = dp sign(p) F(j), x) , 

$M = momentum activity = dp \p\ F(p, x) , /o ^^x 

^ » [ A sign(p)p2 

Qe = energy flux = dp r(j),x) . 

I 2m 



Note that for the stationary Maxwellians of (3.13) these fluxes are equal to 



2rmr 



^p=A/^erf(aV/3m/2), 



$M =—e-^^^ /2 + aW— -erf(aV/3m/2) 



Also note that for a = the quantity <1>m does not vanish. This is because it 
measures the total outgoing flux, not the directed outgoing flux (which is of course 
when a = 0). 

Lemma 3.1. For every stationary solution of (3.5) the 3 fluxes of (3.14) are inde- 
pendent ofx£ [0, 1]. 

Proof. From (3.5a) we deduce that 

dx / dpsign(j))F(j),x) = dpdP (^F(p,x) g(P,x) - F(j),x) g{P,x)j , 

which vanishes since dpdP = dpdP. Similarly, multiplying (3.5a) by p and inte- 
grating over p, we get 

dx / dp \p\ F(p, x)= dpdP p (f(p, x) g(P, x) - F(p, x) g(P, x) 

Multiplying (3.5b) by P and integrating over P, we get 

= dpdP P (f(p, x) g{P, x) - F(j), x) g{P, x) 

Adding these two equations, we see that 

9x / dp IpI F(p, x) = dpdP (p + P) i F(p, x) g{P, x) — F(p, x) g{P, x) 
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But this vanishes, since P + p = P + pby momentum conservation, and using 
again dpdP = dpdP. In a similar way, we first have 

d^ j dp^^F(j>,x) 

dpdP — [F{p,x) g{P ,x) - F{p,x) g{P,x) 

Finally, multiplying this time (3.5b) by P"^ /M, integrating over P, adding to the 
above equation and using energy conservation, we get 

d,jdp\^F(p,x) 

'^^'^^ (ik + ^j {F^^^^9iP,x)-F(p,x)giP,x) 
0. 



Thus, all three fluxes are independent of x, as asserted. D 

4 Formulating the heat-conduction problem 

Based on the stationarity equation (3.5), we now formulate the problem of heat 
conduction in mathematical terms. We imagine that the input of the problem is 
given by prescribing the incoming fluxes on both sides of the system. The system, 
in its stationary state, should then adapt all the other quantities, F and g, to this 
given input, which describes really the forcing of the system. In particular, the 
distribution of the outgoing fluxes will be entirely determined by the incoming 
fluxes. 

We now formulate this question in mathematical terms: The incoming fluxes 
are described by two functions ^o(p) (defined for p > 0) and ^i(p) (defined for 
p < 0). These are the incoming distributions on the left end (index 0) and the right 
end (index 1) of the system. 

In terms of these 2 functions, the problem of existence of a stationary state can 
be formulated as (recall that the rescaled system has length one): 

Is there a solution (F, g) of the equations (3.5) with the boundary conditions 

F(p, 0) = ^o(p) , Vp > and F(p, 1) = ^i(p) , Vp < . (4.1) 

Assume for a moment that, instead of the boundary conditions (4.1) we were 
given just F(p, 0), but now for all p G M, not only for p > 0. Assume furthermore, 
that g(j), x) is determined by (3.5b). In that case, the relation (3.5) can be written 
as a dynamical system in the variable x: 

d,F{-,x) = ^{F{-,x)). (4.2) 
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Thus, if F(-,0) is given, then, in principle, F{-, 1) is determined (uniquely) by 
Eq.(4.2), provided such a solution exists. We denote this map by ^: 

%:F{;0)^F{;1). 

What is of interest for our problem is the restriction of the image of ^to functions 
of negative p only, since that corresponds to the incoming particles from the right 
side, and so we define 

^{F{-, 0))] ip) = e^ip) ■ (%{F{-,0))] ip) = 9^{p) ■ F{p, 1) . 



Using this map ^, we will show that when F(-, 0) varies in a small neighbor- 
hood the map ^ is invertible on its image. By taking inverses the problem of heat 
conduction for our model will be solved for small temperature and flux difference. 

Of course, this needs a careful study of the function space on which ^ is 
supposed to act. This will be done below. 

To formulate the problem more precisely, we change notation, and let 

F+(p) = e+(p)F{p,0) , 
F^(p) = e-{p)F{p,0) , 

F+(p) = e+{p)F{p, 1) , 
Ff (p) = e.{p)F{p, 1) . 

We assume now that Fq^ is fixed once and for all and omit it from the notation. 
Then, we see that ^ can be interpreted as a map which maps the function Fq~ to 
F^ , and we call this map <1>. 

We will show below that for Fq in a small neighborhood V the map $ is 1-1 
onto its image $(2?) and can therefore be inverted. For any F in ^{V), we can take 
Fq = <I>^^(F), and we will have solved the problem of existence of heat flux. 

5 The g equation 

We start here with the study of existence of g for given F. Since (3.5b) does not 
couple different x, we can fix x. The equation (3.5b) is then equivalent to 

g = ^Fig) , 

where the operator Af acting on the function h is defined by 

/ dp F{p{p, P), x)h{P{p, P), x) 



AF{h){P,x) 



J dpF{p,x) 



(provided the denominator does not vanish). Note that for fixed P and p, we can 
solve the collision system (3.2) to find the corresponding P and p, namely 

P = - P p and p = P p . 

Q Q Q Q 
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The action of the operator Af can then be rewritten as 

A ^J,v6 ^ g~^ JdpF(p,x)h{P/Q-{l + q)p/q,x) 

AF(h)(P,x) = „ . (5.1) 

J as F{s, x) 

In order to study this operator notice that it does not depend explicitly on x. It 
is convenient to study instead a family of operators indexed by functions if of the 
momentum only. We define (assuming the integral of (p does not vanish) 

Qj dp ip(p) 
A final change of variables will be useful when we study ^^: 

6 The mathematical setup and the main result 

Having formulated the problem of existence of the stationary solution in general, 
we now fix the mathematical framework in which we can prove this existence. 
This framework, while quite general, depends nevertheless on a certain number of 
technical assumptions which we formulate now. 

We fix once and for all the ratio /i = m/M of the masses, and assume, for 
definiteness, that n G (0, 1). It seems that this condition is not really necessary, 
and probably the condition m ^ M (and the masses non-zero) should work as 
well, but we have not pursued this. 

We next describe a condition on the incoming distribution, called F in the 
earlier sections. The basic idea, inspired from the equilibrium calculations, is that 
F(p, x) should be close to 

i^reference(p) = exp(-/3p^/(2m)) = exp(-ap^) , 

while the derived quantity g(p, x) should be close to 

5reference(p) = exp(-/3p^/(2M)) = exp(-;Uap^) . 

Upon rescaling p, we may assume henceforth that a = 1. 

The operators of the earlier sections will now be described in spaces with 
weights 

y\^u{p) = expi-up"^) , 

where we will choose u = lior the F and u = fiior the g. 
We recall that the operator ^p in "flat" space is 

JdqF{£^)g{q) 



{l + Q)jdqF{q) ■ 
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We then define the integral kernel in the space with weights exp(— p^) for F and 
exp(— /up^) for g, and write 

2 2 

F(p) = e~^ -uCp) , g(j)) = e~^^ n(p) . 

Here, ^ = m/M = (1 — q)/{1 + q), as before. Expressed with u and i; the operator 
^F takes the form 

^Kvu){p) = \ -^— • {L,u)(p) , (6.1) 

{1 + 0) J dqe-1 v{q) 



where 



and 



(L^u)(p) 



Jdqv(a)K(p,qMq), 



K(p, q) = >Vi(M) • >V^(g)/>V» . 
A simple calculation shows that 

K(p, q) = e-(^P-'?)'/(i+^)' . (6.2) 

Our task will be to understand under which conditions the linear operator ^p has 
an eigenvalue 1 . This will be done by showing that K^ is quasi-compact. It is here 
that we were not able to give reasonable bounds on K{p, q) in the case of different 
exponentials for p > and p < 0, which represents different temperatures for 
ingoing and outgoing particles. 

6.1 Function spaces 

We now define spaces which are adapted to the simultaneous requirement of func- 
tions being close to a Gaussian near \p\ = cxd and u and v having limits, and K^ 
being quasi-compact. We define a space ^i of functions u with norm 

Iklki = / e'^''P'\u(p)\Ap. 
Similarly, ^i is the space of functions v with norm 

11^11.^1 = / ^~^ \v{p)\dp ■ 

Thus, the only difference is the absence of the factor ^ = (1 — g)/{l + g) in the 
exponent. 

We also define a smaller space W2, contained in Wi, with the norm 



\du(p)\ + / e~^'P^\u(p)\dp , 
and the analogous space ^2 contained in ^1 with the norm 

\dv{p)\ + / e^^ \v{p)\dp , 



\V\\^2 
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Remark. To simplify notation we write J |dn(p)| instead of the variation norm. 
However, the "integration by parts" formula would hold with the "correct" defini- 
tion of variation as well. 

Lemma 6.1. One has the inclusion ^2 C L°°, and more precisely 

||^i||L°° < / \M +e^ e~^'P'\u(j))\dp < e^'Wuy^ . 

Furthermore, ifu& ^2. then limp_^-|-oo u(p) exists. The maps u h^ linip^-i-oo u(p) 
and u ^^ f dp exp(— ^up^) • u(j)) are continuous functions from ^2 to M. The unit 
ball of'^2 is compact in ^1. 

Analogous statements hold for the spaces ,^2 (defined without the factor n). 

Proof. The first statement is easy, but it will be convenient to have the explicit 
estimates. We have 

u{y) — u{x) = du , 

and therefore 

\u{x)\ < / |dn| + / \u{y)\dy < / \du\ + e^ / e~^''^\u(p)\dp . 

The second statement follows at once since the functions in ^2 are of bounded 
variation. 

For the last assertions, it follows from the inclusion in L°° that the unit ball 
of ^2 is equi-integrable at infinity in L^(e^^*' dp). Moreover, a set of uniformly 
bounded functions of uniformly bounded variation is compact in any h^{K, dp) for 
any compact subset K of M (see [1], Helly's selection principle). D 

6.2 A cone in ^2 

We will work in the space .^2 but we will need a cone (of positive functions, with 
adequate decay) in this space, in order to prove quasi-compactness of K^. 
We define a cone Cj^ in ^2 by the condition 

C,^ = \v £^2, w > and Z • lim v(j))< l\ , (6.3) 

I p^itoo J 

where 

Z = Z{v) = ^ . (6.4) 

J dp e P v(p) 

Lemma 6.2. The cone Cjr has non empty interior (in ^2) cmd is convex. 

Proof. By Lemma 6.1 the maps v 1-^ limp_>-|-oo v{p) and v ^^ f exp(— p^)t;(p)dj) 
are continuous in ^2 and hence the assertion follows. D 

Remark. Note that a function in the interior of the cone is necessarily bounded 
away from zero, since at infinity it must have a non-zero limit and in any compact 
set, if it is never zero, it is bounded away from zero. 

Remark. Note that the function v = I (the Gaussian) is not in the cone C,^. In 
fact, we require that linip^-i-oo F{p)e^ ■ f e^P dp' / J F{p')dp' < 1. 
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6.3 The main result 

On the set C,^, we consider now the spatial evolution equations (3.5) in the vari- 
ables V and Uy (which is the solution of K^u = u with K^ defined in (6.1)): 



sign(p) / dP ( {Wi-v)(p,x) (>V^-u.(,.))(P,x) 



(Wrz;)(p,x)(>V^-u.(,.))(P,x), 

(6.5) 

= sign(p) / AP{Wvv){ -Qp + {1- q)P, x) {Wf,-u^(.,^)){{l + q)p + gP, x) 
- signip) v(p,x) / dPe~''^ Uv(-,x){P,x) , 

with initial condition v(-,x = 0) G C^. We will give a more explicit variant in 
(7.9). 

Any solution of this equation is a function of p and x, and it is easy to verify 
that it satisfies the equation (3.5a). Together with the definition of n^ we have a 
complete solution of the nonlinear system (3.5). Here we assume of course that the 
r.h.s. of the above equation is well defined as a function, so that we can multiply 
by sign(p). 

Theorem 6.3. For any vq G C^, there are a number x^,,, > and a neighborhood 
X)o ofvQ in Cy such that the solution of (6.5) exists for any initial condition vq = 
f (p, 0) G X;o and for any x in the interval [0,Xt,,j]. The function vq h^ x^,,, is 
continuous from Cj^ to M"'' compactified at infinity. We denote by ^x the semi- 
flow integrating (6.5). For any x G [0,Xt,,)], the map ^x '■ v >-^ ^x(v) is a local 
diffeomorphism, i.e., a diffeomorphism on Y^^. 

Note that this implies in particular that the probability densities for g{-,x) and 
F{-,x) remain positive for all x G [0, x^,(,], which is of course crucial from the 
physics point of view. 

We will prove this in Sect. 7 (and in the appendix). 

7 Bound on the operator K^ and proof of Theorem 6.3 

These bounds are the crux of the matter. They actually show, that, under the con- 
ditions on ^2 and the set Cjr, the operator K^ is quasi-compact. In terms of the 
physical problem, this means that the scatterer is not heating up if the incoming 
fluxes are in ^2- 

The object of study is, for v G C^, the operator 

(^-^)(P^ = n^ ^r\'^ ^ ^^ [ v{a)K(p,q)u{q)dq . 
[1 + 0) J e~9 v{q)dq J '-^^ 

and we are asking for a solution u of the equation K^{u) = u. 
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Lemma 7.1. Ifv>0, K^ is a positive (nonnegative) operator, and 
I dp e~^P' {K^u)(p) = I dp e-^'P\(p) 



and 



\K W — 1 



Proof. Easy, compute and take absolute values. Alternately, consider that the prob- 
ability is conserved in the original space. D 

Our main technical result is 

Proposition 7.2. For v G C,^, there exist a C„ () < C, < 1 and an R > (both 
depend on v continuously) such that for any u G ^2 one has the bound 



j \dK^{u)\ <C I \du\+R\\u\y, . 



Proof. Since v will be fixed throughout the study of K^, it will be useful to intro- 
duce the abbreviation Q = Qv for the normalizing factor 

Q = ? — 2 • (7-1) 

(1 + q) f e-'i\{q)dq 

We will use a family of smooth cut-off functions xl (L > 1) which are equal 
to 1 on [— L + ^,L— ^] and which vanish on |g| > L + ^. Let B be a (7°° function 
satisfying < 9 < 1, with Q{q) = for g < -i and Q{q) = 1 for g > i. We 
define xl by 



Xdq) 



'Q{q + L) iiq<-L+l, 
1 ii-L + ^<q<L~^, 

[e{L-q) iiq>L-^. 



The functions xl are C°°, satisfy < xl ^ 1 and ||XlIIl°° ^^ independent of 
L. Let Li and L2 be two positive numbers to be chosen large enough later on 
(depending on v). We will use the partition of unity 

i = XL + xi- 

Using this decomposition of unity with L = Li and L = L2, we write K^ = 

K(,i) + KP + KP with 



(i^«n)(p) = Qj dqv{P-^)K{p,q)u{q) ■ xlM) , 

«2)n)(^) = Q j dqv{P-^)K(jp,q)u{q) ■ xiS^)XL,(p - Qq) , 

{Kf\)(p) = Q [ dqv{I^)K(j>,qMq) • xiM^xi,(j> " QQ) ■ 
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We will now estimate the variation of the three operators separately. 
For the variation of the first term, we find 

d{K^^^u)(p) = + Qy"dgd^(M)_l_ . K(j,,q)u{q)xLM) 
+ QdpJdqv{I^)dpK(p,q) ■ u{q)xLM) ■ 

Using the explicit form of Kip, q) (see Eq.(6.2)), and some easy bounds which we 
defer to the Appendix, we get the bound 



L +- 
f\d{Ki'\){p)\ < 0{1)Q (\\v\\i^^ + f \dv\] j ' \ \u{q)\dq 



1 



(7.2) 



< 0{l)Q { ||w||l- + / |dt;| J e'^^^i+a) / e-^"^ \u{q)\dq 

< const. ||n||5fj • ||f 11,^2 • 

The variation of K^^^ leads to three terms: 
d(Kfn)(p) 

= Qdp dqv{ 2^ ) K(p, q)u{q) ■ Xli (q) x'l2 (P - QQ) 

+ Q j dqdv{P-^)-^^ ■ Kip,qMq) • xUq) XL.ip - Qq) 

+Q dp / dq v{^) dpK(p, q) ■ u{q) ■ xiM^I ^L^^JP - QO) 
■ = d J21 + d J22 + d J23 . 
In these terms, the variables p and q are in the domain 

^ = {(P, g) G IK^ : \p- Qq\ <L2 + \ and \q\ > Li - \] , 
and for Li = 3^L2/(1 — Q^) and L2 sufficiently large we have from Lemma A.3: 
Kip, q) < exp {-Ciigp - qf - C2LI) . (7.3) 

Therefore, we get for dJ2i: 

/|dJ2i| <const.e"^2Li||^||^||„||^ f dpdqe-^^^^^-'i^\ (7.4) 

The integral exists and is uniformly bounded in L2 (since \qp — q\ -^ 00 when 
\q\ -^ CO). 

The term d J22 is handled in a similar way and leads to the bound 



/ 



|dJ22| < const.e"'^2i2||u||oo / \dv\ . (7.5) 
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The following identity is useful: 

dpK(p, q) = -gdqKip, q) . (7.6) 

For the term dJ23 we observe that from (7.6) one gets, upon integrating by 
parts, with the notation 



dJ23 =gQdpJdqv{S^)Kip,q) • 8, [xiM)XL,(p - Qq)) 

+ gQdp j Av{^-^)-^^K(p,q)u{q) ■ X^M^XL.ip - QO) 

+ gQdp v{^^)K(j), q)du{q) ■ XlSi^XL2(P - Ql) 

'■= d J231 + dJ232 + d J233 • 
All these terms are localized in the domain V. In d J231 there appears a derivative 



X =dq [XLSQ)XL2iP - Qq) 
= - Xiiil^XLiiP - qq) 

- qxlM^x'li^p- eq) 

■.=Xi + X2 . 

The terms involving Xi and X2 can be bounded as d J21 and d J22 by observing 
that suppx'j^ C {Ig'l < Li + i}, and similarly for X2. 

The terms d J232 and d J233 are bounded similarly. 

Together, these lead to a bound 

/ \dK^^\u)\ < con?,t.e-'^^^^\v\y^\\u\y^ . (1.1) 

Remark. Note that in this term, the norm \\u\\c^^ appears with a small coefficient, 
while in (7.2) it was \\u\\cg^ (with a large coefficient). 

Finally, we estimate the total variation of K^^\u) and here, the nature of the set 
Cj^ will be important. We have 

d{Kfu){p) 

=Q j dqAv{^-^)^^K(p,q)u{q)-xiS^)xi,(p-Qq) 

-Qdp / dqv{^^^)K(p,q)u{q) ■ XlSi'' Xl^M - QQ) 

+Qdp j Aqv{^-^)dpK(p,q) ■ u{q) ■ xiS^)xi,(P " Qq) 

:=dJ3i +dJ32 +dJ33 . 

The critical term is dJss, but we first deal with the two others which are treated 
similar to earlier cases. 
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For the first term we have by Lemma A.4 which tells us that K is exponentially 
bounded on V: 

/ IdJsil < const. ||n||L°o / IdwCs)] . 

where V is the domain 

V = {(j), q) : \q\ > Li and \p - Qq\ > L2} . 
For the second term, we have, again by Lemma A.4 below, 

|dJ32| < const. IIuIIloo ||f ||l°o 

The last term is more delicate, and uses the property Z ■ linip^-too v(j)) < 1 of the 
definition of the cone vcone, Eq.(6.3). integrate by parts as before using (7.6) and 
get 

dJ33 =gQdp J dv{E^)j:^^K(j),qMq) • xiM^ xi,(j> - QQ) 
gQdp / dqv{^^)K(p, q)u{q) ■ dq{xLi{q)XL.M " 60)) 



qQAp J vi^)K(p, q)du{q) ■ x^q) xi,(P " QQ) 
'■= d J331 + d J332 + d J333 . 

The term dJ33i is bounded like dJ^i. 

In a similar way d J332 and d J32 are bounded by the same methods. 
The term dJ333 makes use of the limit condition in C^. Consider the integral 
of |dJ333|. This leads to a bound and setting L2 = {L2 — ^)/(l + q)'- 



|dJ333(p)| <qQ sup \v{s)\ ■ I \du{q)\ 



■ / dp K(j), q) ■ XlSi^ Xl2(P - 61) 
<gQ I sup / dpK(p,q) ■ sup^ \v{s)\ / \du\ (-yg^ 



Jg|>Li-t- / l^l>^2 



sup \v{s)\ / |du| 



IT 



fdqe 1 v{q) |s|>l^ 

= Z ■ sup \v{s)\ / |dn| , 
s|>L'2 J 

where Z was defined in Eq.(6.4). Collecting all the estimates, we get 

/ |di^,(n)| < C f e~'"^'\u{q)\dq + ((^2) / |dn| 
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where 

C(L2) = 0(1)6-^2^2 11^11,^2 + 0(1) / \dv{s)\ +Z- sup \v(s)\ . 

J\s\>L'^ \s\>L'^ 

Since v belongs to Cj^, it follows that 

lim C(^2)<1, 

L2— »oo 

and the Lemma follows by taking L2 large enough. D 

Proposition 7.3. For any v G C,^, the equation K^{u) = u has a solution in ^2- 
This solution can be chosen positive, it is then unique if we impose \\u\\cg^ = 1. We 
call it u^. The map v ^^ u^ is dijferentiable. 

Proof. We apply the theorem of lonescu-Tulcea and Marinescu [4] to prove the 
existence of u. Since for « > 0, the operator K^ is positivity improving, it follows 
by a well known argument, see e.g., [6] that the peripheral spectrum consists only 
of the simple eigenvalue one and the eigenvector can be chosen positive. If normal- 
ized, it is then unique. Since the operator K^ depends linearly and continuously on 
V (in ^2). the last result follows by analytic perturbation theory (see [5]). D 

We next consider the equation (6.5) for v. 

d^v(p) = sign(j9) (- I e(i-^')(P-9/(i+^»Ve%((i _ ^)^ _ p)/^) u„(g)dg 



-vip) e-^«n,(g)dgl . (7.9) 

Proposition 7.4. The r.h.s. of the equation for v is a C^ vector field on ^2- 
Proof. This follows easily from the fact that the map v ^^ u^\s,C^ . D 

Theorem 7.5. Let vq £ C,^, and assume that vq is bounded below away from zero 
and has nonzero limits at ±00. Then there is a number s = s{v) > such that the 
solution of equation (7.9) with initial condition vq exists in ^2 cmd is nonnegative 
(moreover, it belongs to Cj^j. 

Proof. Follows at once from the previous proposition and the fact that vq is in the 
interior of Cjr. D 

The proof of Theorem 6.3 is now completed by observing that the map $ : 
Vq <-^ ^(vq) is indeed a local diffeomorphism, since it is given as the solution of 
an evolution equation. 



22 8 Remarks and Discussion 

8 Remarks and Discussion 

8.1 The behavior of the solution at p = od 

Consider the limit p ^ cx) in the expression for K^. We need QP — q = 0(1) 
otherwise the Gaussian gives a negligible contribution. In other words, q ~ gp, 
and we are going to assume from now on that ^ > (the other case can be treated 
analogously). This implies p — gq ^-^ (l — g'^)p which also tends to infinity (the 
same infinity). Therefore, 

\/vr f (±cx)) u(±cx)) 

At,n(±oo) = J. 2 • 

J e P v(p)dp 

In particular, if K^u = u and since we assumed 

^/^T V{±(X)) ^ 

f e~P v(p)dp 

we get u{±oo) = 0. 

For the v equation, we have for large p,q ^ p{\ + g) and (1 — g)q —p^ —Q^P- 
Therefore (inverting limit and derivative) we get 



dxv{-^oo) = sign(ibcx)) 



\/tt \ n(±oo) f(=Foo) — f(±oo) I e ^'^ u{q)dq 

V 1-^ J 



Note that the first term vanishes since u(iboo) = 0. Since the integral C(x) 
f e^^'^ u{q, x) is positive, we conclude that formally, 

92:-u(lboo) = ^V{±00)C{X) . 

8.2 Essential spectrum 

Conjecture. The essential spectrum of K^ is the interval [0, a(v)] with 



^/^Tv{±oo) 
a{v) = max ■ 



J e P v(j))dp 

If a(v) < 1 we are looking for an eigenvalue 1 outside the essential spectrum, 
which is the case we have treated. If (t{v) > 1 we would be looking for an eigen- 
value 1 inside the essential spectrum which would be a much more difficult task, 
since it may well not exist. 

Idea of proof: Similar to the above estimates, the operator K-^ should be written as 
something small plus something compact plus something whose essential spectrum 
can be computed. This last part is likely to be the limit operator at infinity. 
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8.3 Dependence on N 

It should be noted that the equation for dxF has, in fact a scahng of the form 

N-^d^F = 0(1) + OiN-^) . 

This means that in the main theorem (Theorem 6.3), the limit Xt,,j of x for which we 
have a result is quite probably bounded by a quantity of the form 1/(A'^ • A(t'o)), 
where A{vq) measures the deviation of the initial condition vq from a Gaussian. 
Thus, either x^g is very small when N is large, or one has to take vq very close to 
a Gaussian. 

Another way to look at this scaling is to introduce a scattering probability 
7 = b/N where 6 > is a constant independent of N. In other words, a particle 
entering the array of cells from the left has for large A^ a probability e~^ to traverse 
all the N cells (and leave on the right) without having experienced any scattering. 
This is analogous to a rarefied gas. It is easy to verify that equation (1.4) is mod- 
ified by a factor b/N multiplying the right hand side, and hence equation (1.5) is 
unchanged. The stationary equations (3.5) become 

5/i:(t,p) = 0(-j5)(i-A^M/(t,^) 



and 



S/Rnt,P) = 0(+p)(l-A)M/(t,^) 



+ TT— ^ / dpg{t-mL/\p\/i^)^ fit- ^^L- ^^L,p) . 



N l-Q 

Equation (3.5a) follows as explained in Section 3 after a rescaling of space by a 
factor b. 

8.4 Discussion 

The model presented in this paper has the nice property that one can control the 
existence of a solution out of equilibrium. In particular, this means that there is no 
heating up of the scatterers in the "chain", when the system is out of equilibrium. 

The reader should note, however, that the initial condition at the boundary, does 
not allow for different temperatures in the strict sense, only for different distribu- 
tions at the ends. For example, a function of the form 



^(P>0) > . , 2 



exp(— ap^), if p > 0, 
exp(— a'p^), if p < 0, 



with a / a' is not covered by Theorem 6.3. The reason for this failure is that we 
could not find an adequate analog of Lemma A.4 for initial conditions of this type, 
and therefore the bounds on the kernel K(p, q) are not good enough. 
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A Appendix: Bounds on K(p, q) 

We study here the kernel K of (6.2), which equals 

with 

E(p, q) = /i/ _ nq^ _ (M)2 = _(^p _ qfj^i + gf . (A.l) 

Lemma A.l. Assume \q\ < L. There are constants C = C{L, g) and D = 
D{L, g) > such that for all p, 

K(j,,q)<Ce-''P\ (A.2) 

and 

\dpK(p,q)\<Ce~''P\ (A.3) 

Proof. Obvious. D 

Lemma A.2. Assume \p — gq\ < L. There are constants C = C{L, g) and D = 
D{L, g) > such that for all q, 

K(j>,q)<Ce~'''^\ (A.4) 

and 

\dpK(p,q)\<Ce-''''\ (A.5) 
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Proof. The proof is as in Lemma A.l, with the difference that now \p — gq\ < 
L. D 

Lemma A.3. Consider the domain D defined by 

^ = {(P, g) G M^ : \p- gq\ <L2 + \ and \q\ > Li - \} , (A.6) 

with 



Li = rr^L2 . (A.7) 



3^ 

For fixed g G (0, 1) and sufficiently large L2 there are positive constants Ci 
and C2 such that for (p,q) & D one has the bound 



K(p, q) < exp {-Ciigp - qf - C2LI) . 

Proof. From the definition of V and (1 — g'^)q = {gp — q) — gip — gq), we find (for 
sufficiently large L2): 

\gp-q\ > {l-g^)\q\-g\p-gq\ > (l-£>')(Li-i)-^(L2 + i) > gL2 . (A.8) 

Using the form 

(QP - qf > \(gp - qf + jg^Ll , (A.9) 

the assertion follows immediately. D 

We next study the region 

^' = {(P, q) ■■ \q\ > Li and \p - gq\ > L2] . (A.IO) 

In this region, we have the obvious bound 
Lemma A. 4. For (p, q) G T>', one has the bound 

E(p,q) = -{^f- 



